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Abstract 

o 

We discuss the role of so-called "intrinsic-charm" operators in semi-leptonic S-meson 
decays, which appear first at order in the heavy quark expansion. We show 

• ^ by explicit calculation that - at scales /i < - the contributions from "intrinsic- 

charm" effects can be absorbed into short-distance coefficient functions multiplying, 
5^ for instance, the Darwin term. Then, the only remnant of "intrinsic charm" are 

logarithms of the form ln(m^/m^), which can be resummed by using renormalization- 
group techniques. As long as the dynamics at the charm-quark scale is perturbative, 
as{mc) <C 1, this implies that no additional non-perturbative matrix elements aside 
from the Darwin and the spin-orbit term have to be introduced at order l/mf. 
Hence, no sources for additional hadronic uncertainties have to be taken into account. 
Similar arguments may be made for higher orders in the 1 / expansion. 
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1 Introduction 



The heavy quark expansion (HQE) has turned out to be a valuable tool for precision 
calculations of heavy hadron decays [1-4] . In particular, due to the HQE for semi-leptonic 
decays, where the 6 — > c transition is described in the framework of a standard local OPE, 
the relative uncertainty in the CKM matrix element \Vcb\ could be reduced to a level below 

2% [5-8]. 

The expansion in inverse powers of the heavy quark mass nib can be set up for both, 
the lepton-energy spectrum as well as for the total decay rate. The non-perturbative 
input, entering the theoretical description, is given by forward matrix elements of local 
operators in the OPE. The leading term represents the partonic rate and does not contain 
any unknown hadronic matrix element. The pcrturbative corrections to the partonic rate 
have been calculated to order a^, recently [9, 10]. Terms of order 1/mf, vanish due to heavy 
quark symmetries. At order I /ml, two hadronic parameters /x^ and Hq appear, which can 
be interpreted as the kinetic energy and the chromomagnetic moment of the heavy quark 
inside the heavy hadron. The short-distance contribution to the coeffcient of /i^ is known 
to order [H]) while that of fiQ is known at tree level. The dimension-6 operators at 
order l/mf define two additional parameters, which correspond to the Darwin term p^, 
and the spin-orbit term p^S' known from the usual non-relativistic reduction of the Dirac 
equation. The coefficients at that order are only known at tree level, so far. The terms at 
order 1/m^ have also been classified, and introduce five new hadronic parameters [12]. 

It has also been pointed out that at order 1/m^ a dimension-6 operator appears, whose 
matrix element could be interpreted as the "intrinsic-charm" content of the B-mcson [13, 
14]. An order-of-magnitude estimate for the effect has been given in [13], and the additional 
uncertainty from the poor knowledge of these matrix elements has been included in the error 
budget for |Vc6| [15]. However, as we are going to show in this paper, the inclusion of an 
"intrinsic-charm" contribution requires a proper definition of the short-distance functions 
appearing in the lepton-energy spectrum, since the "intrinsic-charm" operators and, for 
instance, the Darwin term mix under renormalization. As long as the strong dynamics at 
the charm-mass scale is treated perturbatively, the effect of "intrinsic charm" can entirely 
be absorbed into short-distance coefficients defined at a low hadronic input scale, and 
the non-analytic dependence on the charm-quark mass can be resummed by standard 
renormalization-group techniques, extending the results in [16]. In this case, no additional 
hadronic uncertainty due to "intrinsic charm" has to be included. On the other hand, 
treating the charm-quark as non-perturbative, the hadronic matrix elements of intrinsic- 
charm operators would remain as unknown parameters. In this case, however, the charm- 
quark dependent terms in the standard expressions for the lepton-energy spectrum and the 
total rate have to be modified accordingly, in order to avoid double counting. 

In this paper, we are going to present a systematic study of how "intrinsic-charm" 
effects will enter the theoretical expressions for the lepton-energy spectrum, depending on 
the treatment of the charm-quark mass scale, with particular emphasis on the mixing of 
the "intrinsic-charm" operators into the Darwin term. 
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2 Calculation of the Charm Contribution 



Starting point for the calculation of inclusive rates within the OPE is the hadronic tensor 
W^i, as it appears in the differential rate for b — ^ ciue transitions, 

dV = 16TTGl\Vcb\^W^,L>'''d(l). (1) 

Here (i0 denotes the invariant phase space for the lepton-neutrino pair, and the leptonic 
tensor is given by 

L^'' = 2 {p^.pl + p:p^,^ - g^'^ p, ■ p,^ - le^^'^^ V..p) , (2) 

where e°^^^ = — eoi23 = +1- Using translational invariance, the hadronic tensor may be 
cast into the form 

2MbW^, = j d''xe'^"''^-''> {B{p)\K{x)^,Plc{x) c(0) 7^Pl fe.(0)|5(p)) , (3) 

where Pl = (1 — 75)/2 projects onto left-handed fields, = p'^ /Mb is the velocity of 
the decaying S-meson, and hv{x) denotes the heavy 6-quark field with the phase e"""'''''^ 
factored out. Performing the OPE for this matrix element, the product of the two 6 — c 
currents is matched onto a set of local operators at scales of the order of the 6-quark 
mass rrtfe. Now, as far as the charm-quark mass is concerned, one may take different points 
of view [17]: 

1. One may assume that nib ~ 'rric 3> Aqcd, which means that the short-distance 
matching coefficients and the phase space integrals are functions of the fixed ratio 
p = ml/ ml. In other words, one integrates out (hard) quantum fluctuations with 
virtualities of order m^^ and is left with light-degrees of freedom: light quarks and 
gluons, together with the quasi-static 6-quark field in HQET. In a standard renormal- 
ization scheme like MS, operators with charm fields do not appear at scales /i < mc. 
More precisely, such operators would correspond to quasi-static charm quarks, which 
cannot contribute to the considered matrix elements, {B\bv ■ ■ ■ Cstatic Cstatic • • • K\B) = 
0, because of energy conservation, + 2mc + AE'goft > m^. This is in fact the point 
of view that is usually considered in the precision determination of \Vcb\- 

2. One may consider the power counting mb ^ mc ^ Aqcd, and integrate out hard 
6-quark fluctuations at a different scale than the charm quark. In this case, for the 
first matching at the high scale ph ~ ^b one still has to keep the charm quark 
dynamical, and the corresponding "intrinsic-charm" operators appear in the OPEQ 
The renormalization-group for these operators can be used to scale down to the 

^More precisely, the "intrinsic-charm" operators correspond to local operators for semi-hard charm 
quarks, i.e. quarks with all momentum components of order rric. This is to be distinguished from the 
hard-coUinear (jet) modes for the charm-quark which appear in non-local operators describing the shape- 
function region [18, 19] for h civ. 
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semi-hard scale /i^h ~ rric, where the charm-quark is finally integrated out. As 
before, the "intrinsic-charm" operators then match onto local operators built from 
light fields. Obviously, the main difference compared to case 1 is, that logarithmic 
terms In (mc/mb) can be resummed into short-distance coefficient functions [16], while 
the analytic terms should be expanded in powers of rric/mb ~ ^J-^qcd /"rrib ~ 0.3. 

3. Finally, one may assume that 3> > Aqcd- In this case, one cannot integrate 
out the charm-quark effects perturbatively, and is thus left with genuine intrinsic- 
charm operators, whose hadronic matrix elements have to be defined at a sufficiently 
high scale fio, satisfying rrih > fio ^ ^c- Notice that the matrix elements of the 
intrinsic-charm operators contain the non-analytic dependence on the charm-quark 
mass TTic, and consequently the partonic phase-space integration for the calculation 
of various moments of the differential rate has to be modified accordingly, in order 
to avoid double counting. 

In the following we shall discuss the different cases in turn. 



2.1 nib ~ nic > Aqcd 

As explained above, when integrating out both, the hard 6-quark fluctuations and the 
charm quark, at a common scale /i ~ m;,, we are left with operators built from soft flelds, 
only. Thus the only matrix elements appearing at order l/m\ are the Darwin term p|, and 
the spin-orbit term p\g deflned by (we use the convention of [12], but omit the hat over 
Pd,Pls) 

2MbpI = {B{p)\K {tD^){ivD){iD>')K\B{p)), 

2Mb pis = {B{p)\k {iD^){ivD){iD,){-ia^'') K\B{p)) . (4) 

In the charged-lepton energy spectrum one obtains (among others) a contribution of the 
form 

rfT ^Glml^^^ ^,pI\ ^e{l-y-p) 
dy 

For later use, we have only quoted the most singular term in the limit y = 2Ei/mh 1, 
and p = ml /ml (the full expressions are provided in the appendix). Upon integration 
it yields a logarithmically enhanced contribution to the total rate 

where the ellipses denote the contributions from the sub-leading terms in ^ which are of 
order p In p. Similarly, we identify the leading terms in the moments (see appendix) 

((.-.o)")^^^|K.r^{8(l-.o)"lnp + ...}. (7) 



4 192vr3 l^^^l —y + • • • ^ • 
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Note that for rrih ~ rric the logarithm is actually of order one and represents a regular con- 
tribution to the matching coefficient (and therefore the remaining terms in curly brackets 
enter on the same level). Also, the phase space boundary for y is y < 1 — p which is away 
from ?/ = 1 by an amount of order one. 

A similar logarithmically enhanced term also appears in the partonic rate, 



b 



partonic 1927r^ 

and in the related moment, 



KfeP {l-12p' lnp + ...} 



{l-y) =^|K,r{6pMnp + ...}. (9) 

partonic iy/TT'^ 

In contrast to the Darwin-term contribution, the logarithmic term vanishes in the limit 
p ^ 0. Nevertheless, as has been shown in [16], such "phase-space logs" can be resummed 
into short-distance coefficients, as we are going to discuss in the following. 



2.2 nib > nic > Aqcd 

When we integrate out the b quark first at a scale fih ~ rrib and still keep the charm quark 
dynamical, we have to take into account operators with explicit charm quarks until those 
are finally integrated out at the semi-hard scale fish ~ "^c- In addition to the dimension- 
5 and dimension-6 operators, defining p^,Pg PdiPlsi t\ms finds (at tree level) 
matrix elements of the local "intrinsic-charm" operators 

2Mb = (27r)^ 6\q - m,v) {B{p)\{K I.Pl c) {c^^Pl K)\B{p)) 

+ (27r)4 (^S\q-mbv)] {B{p)\{id^K^,PLc) {c^^PlK)\B{p)) 



dqa 

+ ..., (10) 

which can be interpreted as the probability to find semi-hard (i.e. off-shell) charm quarks 
inside the heavy 5-meson. 

Notice, that the power-counting for the semi-hard charm fields [c] = (mc)^^'^ is now 
different from the ones for soft HQET fields [b^] = A^/^ and therefore it may be convenient 



to use a notation as in [19], where the "intrinsic-charm" operators in the first line of (10) 
are suppressed by A'^ = (mc/mb)^, the ones in the second line by A'^, the kinetic and 
chromomagnetic operators by A^ = (A/m;,)^, and the Darwin and spin-orbit term by A^. 
Due to chiral symmetry, only the A*^ "intrinsic-charm" operators contribute to the partonic 
rate for b ciu, related to the p^ Inp term in ([s]). Additional soft gluon couplings to 
semi-hard charm quarks are further suppressed, and this will give rise to the A^ suppressed 
terms p|, In p in ([6]), descending from the A^ "intrinsic-charm" operators 



Let us consider first the matrix elements of the operator in the first line of (10). They 
may be decomposed in terms of two hadronic parameters, Ti(p) and T2(p), 

iA7rY{Bip)\K-f,PLCCj^PLb,\B{p)) = 2MBin{fx)g^, + T2ip)v^v,) . (11) 
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The contribution to the rate of the matrix element of the local "intrinsic-charm" oper- 
ators is concentrated at small hadronic mass mx and in the endpoint of the lepton energy 
spectrum. Performing the tree-level matching at /i = rrif,, we have 



2ric 



dm\ dy 



5{m\)5{l-y)T 



ic 



and 



dV 



IC 



dy 



5{l-y)T 



IC 



with 



■^ic 



' 247r3 



cb\ 



3Ti(mb) 



(12) 



(13) 



On the other hand, the calculation of the matching coefficients for the contribution of 
and to the total rate now has to be performed in the limit rUc <C nib. Notice, that 

the naive limit p — >^ in ([s]) would give ill-defined expressions. In particular, the integral 

over 

0{l-y) 



dy 



i-y 



would be infrared divergent in the lepton-energy endpoint. As we will see, the new IR 
divergence in the phase-space integration, appearing in the limit p — > 0, is related to 



the UV renormalization of the "intrinsic-charm" operators (11). Defining the hadronic 



parameters Ti^2(p) in the MS scheme, we also have to perform the phase-space integral in 
D = 4 — 2e dimensions. As a result, the contribution of the Darwin term to the total rate 
is regularized by plus-distributions. 



0{i-y) 



- 5(1 - y) In 



mt 



(14) 



which exactly subtracts the effects of semi-hard charm quarks, that would otherwise be 



double-counted when adding (13) to the decay rate. 



The final expression for the combined contributions of the Darwin term and the "intrinsic- 
charm" operators to the lepton-energy spectrum at order 1/ml can be written as 



dy 



G^ml |2/C'pd(?/,/^)p?)(^) , C'ti(|/,/^)Ti(p) 
ycb\ s 5 H 



247r3 



my 



mi 



(15) 



which should be used for < p < m;,. The matching conditions for the short-distance 
coefficient functions - including the limit p — > for the sub-leading terms in ([s]) as given 
in the appendix - are given by 



Cpo{y,mb) = - 



y^i9-5y + 2y^)e{l-y) 
6(1 -y) 



+ Y/iy-^)-^S"iy-l) + Oias) 



CT,{y,mb) = -36{y~l) + 0{as 
CT^{y,mb) = 0{as) . 



(16) 
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Figure 1: Leading diagrams determining the mixing of fom'-quark into two- 
quark operators. 



In appendix |B.l , we derive the leading terms in the anomalous dimension matrix that 
describe the mixing of the "intrinsic-charm" operators {Ti(/i), T2{fi)} into the Darwin term 
PD(/i), see also Fig. [T|b), 



d 




"^^f^ V T2(/i) / 1 I 4/3 



0^ 

2/3 1 } I ri(/x) I . (17) 




Neglecting the 0{as) contributions to the anomalous dimension matrix, we only determine 
the leading-logarithmic terms which are generated by the renormalization-group equation 
for the short-distance coefficients 

1 

Cpoiy^fJ') - Cpoiy^^b) - - In — {CT,{y,mb) -2CT2{y,mb)) . (18) 

Now, integrating out the semi-hard charm quarks at = is equivalent to setting 

Ti(/i<m,) = 0. (19) 



In this case, the expression for the lepton-energy spectrum (15) simplifies to 



dy 



pI+ic 247r3 ' ml ' ^ ' 



and the information on "intrinsic charm" , i.e. the non-analytic dependence on the charm- 
quark mass, has been completely absorbed into the short-distance function Cp^{y,mc). 
This can be made explicit by inserting the leading-order matching conditions (16) for 
CTXy^^b), which results in 

2 

TTl 

Cp^{y,m^) ^ Cp^(y,m6) + ln— I 6{y - 1) . (21) 

^Strictly speaking, these are N~-'^LL. To resum the leading-logarithms of the order a" x In" p, we would 
need the 0{as) mixing of the Darwin term into itself and the complete set of "intrinsic-charm" operators 
into themselves, which goes beyond the scope of this work (see however [16] for a complete leading-log 
analysis for the dimension-7 contributions to the total rate). 
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In this way (20) reproduces the logarithmic term in the lepton-energy moments in ([T]) as 
well as the finite terms (given by the limit p ^ of Eq. (39) in the Appendix). 

Similar considerations can be made for the In p term in the partonic rate. We 
decompose the matrix elements of the operators in the second line of (10) as 

(47r)2 {B{p)\{id^K^,PLc) {c^,PlK)\B{p)) 
= 2Mb (Tsip) g^,v + Ti{p) g^^ + T^{p) g^a + Tq{ij,) v^v^v^ - Tyd^) ie^^apv^) . 

(22) 

(Notice that in unpolarized observables, only the sum T^^p) +T5(/i) appears.) Generalizing 
the results for the total rate in [16] to the lepton-energy spectrum, and concentrating again 
on the leading logarithmic terms, we find 



dT 
dy 



partonic+IC 



Vcb?\ Co{y,iJ,) + pCi{y,p) + p^C2{y,p) 



mi 



+ 



(23) 



with 



Co(y, m,) = (6y2 - %3) ^^(1 -y) + 0{as) , 
Ci(y, m,) = -6i/2 ^^(1 -y)-Q5{y-l) + 0(«,) 

C2(y,"^f,) 



\l2e{l-y)^ 




'Qe{l-y)- 




+ 





and 



- 6 ^(1 - ?/) + 6 (5(y - 1) + 3 5'{y - 1) + 0{a, 

Cr.iy, m,) = -2A6'iy - 1) + 48 5(y - 1) + 0(a,) , 
CT,My,^b) = -245(y - 1) + Oias) , 
CTeiy,mb) = 0{as), 
CTr{y,m,) = 246'{y - 1) + 0{a,) . 



(24) 



(25) 



Again, the "intrinsic-charm" operators mix into the 2-particle operator m^hv^'-^^h 
(see appendix), and consequently, the coefficient C2{y,p) evolves as 

,2 



C2{y,mc) ~ C2{y,mb) 



In 



p- 



mi 



{CT.,{y,mh) - CT^{y,mb) - CT,{y,mb) - CT,{y,mb)) . 



(26) 



Inserting the leading-order matching conditions, one has 

{CT^y^nik) - CT^{y,mb) - CT,{y,mb) - CTr{y,mb)) = 66'{y - 1) - 126{y - 1) , 

(27) 

and one reproduces the logarithmic terms — 12p^lnp in Fpart and 6p^lnp in (1 — ?/)part, 
respectively, see ( 8|9 ). 



2.3 nib > nic ~ Aqcd 

If we consider the dynamics at the charm-quark mass scale to be in the non-perturbative 



regime, we cannot exploit the condition ( 19 ) and are left with the general formula for the 



leptonic-energy spectrum in (15), which should be evaluated at a scale fiQ that satisfies 
T^c -C /io < '^b- Moreover, we have to take seriously the new power counting which implies 
that terms of order now count as (A/m^)^ and should be neglected to the order that we 
are considering, we are thus left with 



dr 

dy 



dy 



partonic 



Co{y,fio) + pCi{y,iJ,o) + 0{p'^ 



Vcbl 



247r3 



+ 



mi 



(28) 
(29) 



together with the contributions to the lepton-energy spectrum from yU^,/i^ and Pls (^^^ 
e.g. [12]), where the limit p — >■ to the considered order (1/m^) is trivial. 

In that order, the genuinely intrinsic-charm contribution comes together with the Dar- 



win term, only. In particular, to leading logarithmic accuracy (18), the contributions to 
the total rate, and the moments {y) and (y^) can be obtained as 



r{3) 
(y) 



247r3 
Glm 



F"''b IT/ |2 



247r3 



X(/io) + 



X(/io) + 



ml 



PoiP'o) 



ml 



pUpo) 



ml 




(30) 
(31) 
(32) 



where we defined the parameter combination 

3Ti(/io) 



+ In 



pI pUpo) 



ml 



mt 



ml 



(33) 



Considering a sizeable value for Ti{fiQ) at small hadronic scales (in contrast to the per- 
turbative situation considered in the previous subsection), and taking into account that 
the pI, contribution in X(po) is formally enhanced by In fi^/ml, we may ignore the (small) 
differences between the individual moments induced by the numbers in square brackets 
in ( [30|31[32[ ), to first approximation. Therefore, even in this genuine intrinsic-charm sce- 
nario, the inclusion of a large non-perturbative intrinsic-charm effect, basically amounts 
to treating the Darwin term pfj for the effective parameter X. In any case, one may con- 
sider the limit mc ~ Aqcd rather academic, and would prefer the scenario with semi-hard 
charm-quarks as in the previous subsection for the discussion of "intrinsic-charm" effects 
in inclusive semi-leptonic B decays. 
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provide the appropriate formulas for the massless 
hmit, relevant to 6 — > uiu decays, after appropriate changes Vet Vub and re-interpretation 
of the intrinsic-charm operators as so-called weak annihilation operators [20,21]. Notice 
that the (local) annihilation operators enter at order 1/ml in the standard OPE, whereas 
their non-local counterparts, necessary to describe the shape-function region, already enter 
at (relative) order A /nib [22-24]. 

3 Conclusion 

We have shown how the "intrinsic-charm" contribution in semi-leptonic 5-meson decays 
is related to the renormalization of sub-leading operators (like m'^byby and the Darwin 
term) appearing in the operator product expansion for the lepton-energy spectrum and 
the total rate. We have distinguished three different cases which correspond to different 
power counting for the charm-quark mass. In the first case, one assumes rrib ~ rric, i.e. the 
charm-quark is already integrated out at the hard scale, set by the large 6-quark mass in the 
OPE. Consequently, all dependence on the charm-quark dynamics is already encoded in the 
matching conditions for the hard coefficient functions, and no "intrinsic-charm" operators 
should be introduced below the hard scale. The only remnant of "intrinsic charm" is the 
non-analytic dependence of the coefficient functions on the ratio p = m'^/ml. 

Another viable scenario treats the charm quark mass as intermediate between the hard 
and the soft scale in the OPE, mb ^ rric 3> Aqcd- In that case, four-quark operators 
including soft fe-quark fields and semi-hard charm quarks have to be included in the OPE. 
At the same time, in order to avoid double counting, the semi-hard region has to be 
subtracted from phase-space integrals by a suitable regularization of the decay spectra in 
the limit rric m^. We have shown by explicit calculation how the mixing between the 
"intrinsic-charm" operators and the Darwin term generates the logarithmically enhanced 
terms entering the OPE at order l/mf. Similarly, extending the results of [16], we could 
reproduce terms of order In p in the partonic rate. After integrating out the charm 
quark at the semi-hard scale, the moments of the lepton-energy spectrum can be entirely 
described in terms of the standard hadronic input parameters, whereas - again - the 
complete charm-quark dependence enters via (eventually renormalization-group improved) 
short-distance coefficients, multiplying, for instance, the Darwin term. 

A somewhat more exotic approach would treat the charm quark as light, i.e. of order 
^QCD- Only in this case genuine intrinsic-charm (i.e. non-perturbative) effects have to be 
taken into account. Still, we have found that on the level of a few lepton-energy moments, 
the experimental data basically constrain a particular combination of the intrinsic-charm 
contribution and the Darwin term, such that to order 1/m^ the number of independent 
hadronic parameters effectively remains the same. 

The main conclusion to be drawn is that, as long as the strong dynamics at the charm- 
quark scale can be treated perturbatively, "intrinsic-charm" effects do not induce an addi- 
tional source of hadronic uncertainties at the level of 1/ml power corrections, apart from 
the usually considered Darwin and spin-orbit terms. The same will be true for higher or- 



We should also mention that (28p9 
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ders in the l/rnt expansion as classified in [12]. The issue of whether to resum logarithms 
ln(m^/m^) by introducing the above 2-step matching procedure, or sticking to the stan- 
dard 1-step matching has to be decided by considering radiative corrections to the 1/m^ 
expressions which is beyond the scope of this work (see also [10]). 
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A Lepton-Energy Spectrum and Moments 
A.l Partonic rate 

The complete expression for the partonic contribution to the lepton-energy spectrum with 
m6 ~ m-c > /i is given by [12] 



Wcbl 



4p3 



6 (p3 + p2) I2p2 



partonic 1927r3 ' ' [V (?/ - 1)^ {v - 1)^ y-1 

4y' - 6(p - l)y' + 2(p - 3)p')e{l - y - p) } . 



dT 

dy 



From this one can obtain closed expressions for the (1 — y)" moments, 

12(p"-l)p2 4(p"-p2)p 12(p"+2-l)p 



(34) 



partonic 1927r'^ ^ U — 2 U + 2 

6 (p2 + p) (p" - p) 2 (p3 - 3p2 - 3p + 1) (p"+i - 1) 



+ 
+ 



n — 1 n + 1 

6(p + 1) (p"+3 _ 1) 4(p'^+^-l) 



n + 3 

Arbitrary moments can be derived via 



n + 4 



(35) 



k=0 ^ ^ 



(36) 



Expanding (35) in the small parameter p, the logarithmically enhanced terms at order 



p appear only in the total rate and the first moment 



((i-i/D 



partonic 



" 1927r3 



(12 5"° - 6 5"^) p^ Inp + analytic/h iffher-order terms in 



(37) 
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A. 2 Darwin-term contribution 

The full contribution related to the Darwin term in the lepton-energy spectrum for the 
case rrih ~ > is given by [12] 









IT/ |2 Pd 

\Vcb\ ^ 

40p3 


dy 




1927r3 






( 


3(?/ - 1) 








28p 



^ 8p2(3p+l) ^ 6p2(3p + l) ^ 16p(2p2 - p - 1) 



{y - 1)^ 



3(y-l)3 



3(,-l)2 .■^ + 3(5p^-5p^ + 10p + 22) 
+ ^(p + 3)(y - 1) + 4(y - If + ^{y - lA 9(1- y- p) 



5{l-y-p)\. 



2{p-l)\p + lf 
3p2 

From this one can obtain closed expressions for the moments, 
Glml ,2p1,/8(p"-1) 2 

V r 



((i-i/D 



Pi, 



1927r3 



+ 



^1^^ 3^^-') ^P^^^P + 3(n-l) 

2 (5p3 - 5p2 + lOp + 22) (p"+i - 1) 8(p + 3) (p"+2 _ i) 

3(n + 1) ^ 3(n + 2) 

4(p"+3_i) 8(p"+4-l) 16 (2p2 - p - 1) (p2 - p") 

n + 3 ^ 3(n + 4) 3(ra - 2)p 

6(3p + 1) (p3 - p") 8(3p + 1) (p^ - p") 40 (p5-p") 



[n - 3)p 



{n - 4)p2 



3(n - 5)p2 



(3^ 



(39) 



Taking the limit p ^ in (39), the logarithmically enhanced terms appear only in the 
total rate 



Kfer^(5"° lnp + 0(plnp)) 



(40) 



B Operator Mixing 
B.l Dimension-6 

In the following we briefly sketch the derivation of the elements of the anomalous- dimension 
matrix that govern the mixing of the four-quark ("intrinsic-charm") operators into the 
Darwin term. For simplicity, we do not construct the complete set of independent operators 
that would be needed to describe the full one-loop anomalous-dimension matrix, but rather 
focus on the effect of the charm-loop diagram in Fig.fTtb). For this purpose it is sufficient 
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to consider the two operator structures which enter the hadronic tensor at tree- level (10): 
{Any 



2MbTM = ^^(A{B{p)\KiPLC-cfPLK\B{p)) - {B{p)%^,Plc-c^^PMB{p)) 



3 

(4vr)2 

(41) 

Together with the Darwin term they are used to define a simplified operator basis 
Op^=K {iD^){ivD){iD^^)K, 

Ot, = (Anf fx'' ^ {K I^Pl c c^^Pl K - k ^Pl c c^Pl K) , 

Ot, = {Aiif fi^' ^ (4 K fPL c cfPL K - K I^Pl c cYPl K) • (42) 

Notice that for convenience, we have extracted a factor (47r)^ fi^'^, in order to have a simple, 
dimensionless anomalous-dimension matrix|3 

Calculating the one-loop matrix elements of the operators Ot^ 2 for the partonic tran- 
sition b ^ b in the presence of a soft background field A^{k), see Fig. [l](b), and comparing 
with the tree-level matrix element of the Darwin-term operator, we obtain the following 
results in D = A — 2e dimensions, 

= -I (J (&|(^pJ&)tree, (43) 

where the one-gluon matrix element of the Darwin-term operator on parton level is given 
by 

{b\Op^\b)tree = \ {b\k [tD^, [{iV ■ D),lD'']] b,\b)tree + 0{l/mb) 

= I ((t; ■ k){k .A)-e{v A)) UkUt, + ... (44) 



From (43) we read off the desired elements of the anomalous dimension matrix 

/ 0\ 

7 = -2/3 + 0{as) , (45) 
\ 4/3 0/ 

where the neglected higher-order terms describe the mixing of "intrinsic-charm" operators 
into themselves and of the Darwin term into itself, which are not explicitly needed for the 
discussion in the body of the text. 



■^With this convention, the anomalous-dimension matrix is of order (a^)''. In order to have it in the 
standard form, one would have to extract a factor — Anag /i^*^, instead. 
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B.2 Dimension- 7 

A similarly simplified analysis can be performed for the mixing of the dimension-7 "intrinsic- 
charm" operators into the dimension-7 two-quark operator m^h^fh^. As before, defining 

Ot, = (47r)2 /i^^ i {{tv ■ dk lf.PL c) {cYPl K) - {tv ■ dk fPL c) {c^Pl k)) , 

Ot, = (47r)2 i {{td^ k fPL c) (c7"^L h,) - {%v ■ dk ^Pl c) {c^jPl K)) , 

On = (47r)2 /i^^ i {{td^ k ^Pl c) {cfPL k) - {iv ■ dk ^Pl c) {c^Pl K)) , 

On = (47r)2 /i^^ i (6 {tv ■ dk i^PL c) {ci/jPl k) - (u; ■ dk I^^Pl c) {c^Pl k)) 

- (47r)2 i {{id^ k fPL c) {c^Pl k) + k rPL c) {cfPL k)) , 

On = W fi^' ^ e'^'^"^ k IuPl c) {ci,Pl k) , (46) 

we calculate the contributions to the 2-parton matrix elements from the tadpole diagram 
in Fig. [l]^a) as 



8 \e 



\O2\k 



tree ; 



tree ; 



{b\On\k^'^ = {b\On\k^'^ = f- + In 4 + • • •) {b\0,\k 
{b\On\b)^'^ = , 

{b\On\by'^ = f- + In 4 + • • •) (&|02|&)tree , (47) 

8 \e J 



from which we read off the elements of the anomalous dimension matrix entering (26) 
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